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We have studied nonlinear superconducting resonators: A/2 coplanar- waveguide (CPW) res- 
onators with Josephson junctions (JJs) placed in the middle and A/4 CPW resonators terminated 
by JJs, which can be used for the qubit readout as "bifurcation amplifiers." The nonlinearity of 
the resonators arises from the Josephson junctions, and because of the nonlinearity, the resonators 
with appropriate parameters are expected to show a hysteretic response to the frequency sweep, 
or "bifurcation," when they are driven with a sufficiently large power. We designed and fabricated 
resonators whose resonant frequencies were around 10 GHz. We characterized the resonators at low 
temperatures, T < 0.05 K, and confirmed that they indeed exhibited hysteresis. The sizes of the 
hysteresis, however, are sometimes considerably smaller than the predictions based on the loaded 
quality factor in the weak drive regime. When the discrepancy appears, it is mostly explained by 
taking into account the internal loss, which often increases in our resonators with increasing drive 
power in the relevant power range. As a possible origin of the power-dependent loss, the quasiparticle 
channel of conductance of the JJs is discussed. 
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I. INTRODUCTION 

In many experiments related to quantum information 
processing, the readout plays a critical role. One of the 
recent major steps forward in the experiments of super- 
conducting qubits was related to the improvement of 
readout, and the quantum nondemolition readout was 
reported/ In the readout scheme of Ref. [l], the charge 
qubit was non-resonantly coupled to a superconducting 
linear resonator. The qubit state was detected as a shift 
in the resonant frequency of the resonator. The measure- 
ments were done with a weak driving power of n ~ 1, 
where n is the number of measurement photons popu- 
lated in the resonator. 

It has been suggested that employing a nonlinear res- 
onator instead of linear resonator should relax the strong 
demand of low-noise broadband microwave measure- 
ments because the latching effect between the bistable 
states in a nonlinear resonator is expected to provide a 
larger but still fast enough response. 2 Such scheme, or the 
"bifurcation amplifier," has been applied to the readout 
of charge qubits^^ and flux qubit Sj^i 7 -^ where Joseph- 
son junctions (JJs) as nonlinear inductors are employed 
for the nonlinear resonators. In many cases, the nonlin- 
ear resonators with JJs have lumped-element capacitors, 
and their resonant frequencies are around 1 GHz. 3 >2&&& 
In general, however, using a distributed element such as 
a coplanar waveguide (CPW)^ ' 111 ' 11 ' 12 makes it easier to 
increase the resonant frequency and the quality factor. 
A higher resonant frequency is advantageous for sup- 
pressing photon- number fluctuations, which may cause 
qubit dephasing even when the readout circuit is turned 
offji^ Choosing a high quality factor and a resonant fre- 
quency much higher than the qubit frequency would be 
also helpful in suppressing the qubit relaxation through 



the resonator. 14 

In this work, we have characterized at low tempera- 
tures, T < 0.05 K, nonlinear CPW resonators with JJs, 
whose resonant frequencies are as high as 10 — 11 GHz. 
We observed the "bifurcation" as expected theoretically; 
the resonators showed a hysteretic response to the fre- 
quency sweep when they were driven with a sufficiently 
large power. We analyze the hysteresis in detail, paying 
close attention to the internal loss, which has not nec- 
essarily been examined well in the earlier works&i - on 
superconducting nonlinear resonators designed for read- 
ing out flux and charge qubits. 



II. THEORY 

The nonlinear resonators studied in this work are 
shown schematically in Figs.[IJa) and[TJb). They are A/2 
and A/4 CPW resonators, respectively. The crosses in the 
figures denote either single JJs or dc SQUIDs, the origin 
of the nonlinearity. We treat dc SQUIDs as single JJs 
with tunable critical currents Iq. Here, we should men- 
tion that superconducting transmission lines can also ex- 
hibit nonlinear ityi 1 5 ' 1 6 > 1 7 The microwave power required 
for this effect, however, is many orders of magnitude 
larger than the power range of our experiments. Thus, we 
treat superconducting CPWs as linear elements. Port 1 
in Figs.[lja) andQJb) is connected to a microwave source, 
whose equivalent circuits are Figs.[ljc) and[ljd). For uni- 
formity, we always use Fig. [ljc) in this paper. The A/2 
resonator in Fig. [Ua) nas been mappe d 10 ' 11 to a series 
circuit of LCR and JJ [see Fig. [IJe)] in the vicinity of 
the resonance, by considering the electromagnetic envi- 
ronment seen from the JJ in the linear regime, where 
the drive is so weak that the resonator does not manifest 
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FIG. 1: Schematic diagrams of microwave circuits. (a) 
A/2 and (b) A/4 coplanar-waveguide (CPW) resonators with 
Josephson junctions (JJs), where To and I' are the critical cur- 
rents of the junctions, I and V are the CPW length, and C ln , 
Cout, and C' c are the coupling capacitance to the microwave 
line. Both in (a) and (b), Port 1 is connected to a microwave 
source, which can be expressed as (c) an ideal voltage source 
and Zo — 50 Q in series, or (d) an ideal current source and Zo 
in parallel. Port 2 is terminated by Zo in (a) , whereas shorted 
in (b). In the vicinity of the resonance, (a) has been mapped 
to (e) a series circuit of LCR and JJ . 10 ' 11 



nonlincarity. The series circuit is one of the nonlinear 
systems studied by Manucharyan et al^ They solved 
the equation for the Duffing oscillator by retaining the 
terms with the first harmonic only, then, mapped several 
nonlinear systems to the Duffing oscillator by comparing 
the differential equations of the systems with the Duffing 
equation^ 

The series resonant circuit in Fig.QJe) "bifurcates" and 
becomes bistable in a certain frequency range when the 
drive amplitude Vd exceeds a critical value V c , provided 
that the amplitude I c of the current through JJ at Vd = 
V c is sufficiently smaller than Iq . According to Ref . [l8l 
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is the loaded quality factor, Q e is the external quality 
factor, Q u is the unloaded quality factor, 

L JQ = ipo/Io (4) 

is the Josephson inductance, tpg — %/(2e), u> c satisfies 

n c = 2Q l (l-u; c /u; Q ) = y/3, (5) 

and coo is the resonant angular frequency in the linear 
regime. Equations ^ and suggest that a large Qi 



and a small Iq are favorable for the bifurcation. When 
Ql ^> 1, the normalized boundaries of the bistable region 
are expressed as 
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2(3/(1 — uj/loq) and ui is the angular 



where f2 = 
frequency^ 

The mappin g 10 ' 11 of Fig. QJa) to Fig. QJe) has been 
done for circuits with L 3> Ljq by assuming that 



Vd = ZouJodnV, 



(7) 



where Zq = 50 f2 and C- lri is the coupling capacitance to 
the input microwave line, and that 



L + L Ja ~L = ttZ / (2u/ ). 



(8) 



Note that the mapping is based on the behavior in the 
linear regime. We have simulated the circuit in Fig. [TJa) 
taking into account the drive-strength dependence, as we 
will describe in the next paragraphs, and confirmed that 
the simulated V c and the value calculated from Eqs. |T]), 
, and ((51) agree reasonably well. 

The heart of our simulation is expressing JJ as an 
inductance that depends on the amplitude A of the 
superconducting-phase oscillation, 
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where J is the Bessel function of the first kind. The idea 
is based on the following considerations: suppose that 
the superconducting phase S of JJ is oscillating as 



S = Asinwr;. 



(10) 



Then, the voltage Vj across the junction and the current 
I j through the junction are given by 



Vj/ipo — dS/dt = Aucosut, 



(11) 



and 



Ij/h 



sin 5 = sin(A sin ujt) 

oo 

2^ J 2 fc+i(A)sm[(2fc + rM]. (12) 



fc=0 



Here, let us retain the terms with the first harmonic only, 
as Manucharyan et al^ did for solving the Duffing equa- 
tion. Then, Eq. (fl"2|) reduces to 



(13) 



Ij/I = 2Jx(A) sinLut. 



By comparing Eqs. (|11|) and (|13[) with the equations for 
a usual inductor L, Vl <x toLcosujt when 1^ oc sinwrj, we 
obtain Eq. Note that Lj is nonlinear in the sense 
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FIG. 2: (Color online) Theoretical boundaries of the bistable 
region for a nonlinear resonator. The solid and broken curves 
are the simulation and Eq. (JSJ), respectively. 
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FIG. 3: (Color online) Simulated drive dependence for the 
resonator in Fig. [2] at / = 10.30 GHz. (a) Amplitude of the 
superconducting-phase oscillation [see Eq. () 10[l ] vs. drive volt- 
age amplitude, (b) Phase of the transmission coefficient S21 
vs. incident power to the resonator. The markers in both (a) 
and (b) correspond to the boundaries in Fig. [3] The arrows in 
(b) indicate the expected switching in actual measurements. 



that it depends on A. Its impedance Zj, on the other 
hand, is written as 

Zj=JljLj, (14) 

as if it were a usual linear inductor because we retained 
the terms with the first harmonic only. 

Equations (J9j) and (fT4|) make the simulation easy. Now, 
for example, it is straightforward to include in the calcu- 
lations, if necessary, the junction capacitance Cj and the 
quasiparticle resistance R qp , both of which are parallel to 
Lj. The relationship between 7? qp and the quality fac- 
tors will be discussed in Sec. El We used the transmission 
(ABCD) matrix (see, for example, Sec. 5.5 of Ref. 
for the simulation. The matrix Tz for an impedance Z 
is 

r* = (j?), (15) 

and that for a section of lossless CPW with a length I is 

_/ cos (31 j-ZcpwSin/jA , . 

Jcpw- ^( Zcpw )-i sinjSJ cos/3Z )> [W) 

where j is the imaginary unit, Z cpw is the characteristic 
impedance, f3 = uj/v p , and v p is the phase velocity. When 
CPW has a loss, (3 in Eq. (fTH)) is replaced by (a + j(3)/j, 
where a characterizes the loss. 

An example of the simulation is shown in Fig. [21 It is 
for a resonator of the Fig.QJa) type. The parameters are 
I = 5 /iA, I = 2.8 mm, and C in = C ou t = C c = 5 fF. 
Both CPWs are lossless with Z cpw — 50 il and v p = 0.4c, 
where c is the speed of light. When we neglected Cj 
and assumed that R^s — 0, the simulation yielded 
fa = lu /(2tt) = 10.306 GHz, Qi = Q e = 3.0 x 10 3 , 
and V c = 2.2 /iV. This simulated V c and the value cal- 
culated from Eqs. (P), (J7J), and © agree within 2%. In 
Fig. [2l the solid curves are the simulated boundaries of 
the bistable region, and compared with Eq. ([S]), the bro- 
ken curves. The vertical axis of Figgis the drive voltage 
amplitude V normalized by V c . In general, the solid and 



broken curves are similar. As the drive amplitude is in- 
creased, the difference becomes visible especially for the 
lower boundary, but this trend should be reasonable be- 
cause the mapping of Eqs. J7]) and {SJ is based on the 
behavior in the linear regime as we mentioned earlier, 
and because at a given frequency, A at the lower bound- 
ary is actually larger than A at the upper boundary as 
we will see in Fig. [3ja). 

The simulated boundaries in Fig. [5] were obtained by 
calculating A vs. V at different frequencies. Such calcula- 
tion for / = 10.30 GHz is shown in Fig.[3](a); the marked 
local maximum and minimum of V are the boundaries. 
In experiments, we measure instead of A, the transmis- 
sion coefficient S21 of A/2 resonators and the reflection 
coefficient T of A/4 resonators. These quantities are also 
simulated by the same method. The phase of S21 is com- 
puted in Fig. [3Jb) as a function of the incident power 
Pi n , which is more relevant than V from the experimen- 
tal point of view. The power is given by 

P in = V?/(2Z ), (17) 

where 

= V/2 (18) 

is the amplitude of the incident voltage. Because the 
section between the markers is unstable, the resonator 
response in actual measurements is expected to be hys- 
teretic, that is, switching indicated by the arrows in 
Fig. (2Kb) should be observed. 

Regarding A/4 resonators, when we chose in Fig.fijb), 
Iq = 2Iq, C' c = C c , I' — I, and the same unit-length 
CPW properties, the simulation yielded essentially the 
same ojq, Qi, V c , etc. as those of the A/ 2 resonator in 
Figs, rj and E3 

When we compare the theoretical predictions in this 
section with experimental results, we should note that in 
experiments, we cannot make the drive amplitude arbi- 
trarily large in order to see the bistability arising from 
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TABLE I: List of resonators. A1-A4 have single Josephson 
junctions, whereas B1-B2 have dc SQUIDs. fo is the resonant 
frequency in the linear regime, Qi is the loaded quality factor, 
and P c is the critical incident power for the bifurcation. All 
values are for the zero magnetic field, and Qi was calculated 
in the linear regime. 



Reso- 




fo 


Qi 


Pc 


nator 


Type 


(GHz) 


(xlO 3 ) 


(dBm) 


Al 


A/2 


9.86 


6 


-97 ± 2 


A2 


A/4 


9.70 


3 


-100 ± 1 


A3 


A/2 


9.98 


1.1 


-101 ± 2 


A4 


A/4 


9.83 


1.6 


-103 ± 3 


Bl 


A/2 


11.30 


1.0 


-85 ± 1 


B2 


A/4 


11.12 


1.3 


-82 ± 1 



the nonlinear inductance of JJs because for example, Ij 
has to be smaller than Iq. Let us also recall that both 
Eq. © and the simulation are based on the approxima- 
tion of retaining the terms with the first harmonic only; 
the accuracy of the approximation becomes worse as the 
drive amplitude increases. 

III. EXPERIMENT 

We studied two series of CPW resonators listed in Ta- 
ble HI The fabrication of Series A, which has single JJs, 
was almost the same as that of the sample in Ref. [T3; a 
Si wafer covered by a layer of thermally oxidized Si02 
was used, Nb interdigital coupling capacitors and Nb 
CPWs were patterned by photolithography and reactive 
ion etching, and Al JJs by electron-beam lithography and 
shadow evaporation. In order to realize a superconduct- 
ing contact between Nb and Al, the surface of Nb was 
cleaned by Ar + milling before the Al evaporation. In Se- 
ries B with dc SQUIDs, on the other hand, all electrodes 
including those of SQUIDs are Nb, and everything was 
fabricated by the photolithographic technology. Thus, 
between the two series of resonators, there are a couple 
of differences regarding the CPW that we should note: 
the thickness of Nb film and the quality of SiC>2. 

The Nb film is much thinner in Series A (0.05 /im) than 
in Series B (0.4 /xm). As a result, highly likely due to the 
kinetic inductance,— we saw a noticeable difference in 
v p . In Series A, Nb is deposited directly on the ther- 
mally oxidized Si02 of the Si wafer, whereas in Series B, 
there is a layer of sputtered Si0 2 between Nb and ther- 
mally oxidized SiC>2. It seems that this sputtered SiC>2 
decreases Q u . According to Ref. |2l|, the loss tangent of 
Si02 grown by plasma-enhanced chemical vapor deposi- 
tion (PECVD) is much larger than that of thermally ox- 
idized SiC>2. Thus, it is highly likely that sputtered SiC>2 
is also much more lossy than thermally oxidized SiC>2- 
We will come back to this point in Sec. IIVBI 

Our A/2 resonators (Al, A3, and Bl) are designed to 
be symmetric, that is, in Fig. [TJa), C ln = C out = C c 
and the two CPWs are the same. For all resonators, we 



intended to obtain /o ~ 10 GHz and Qi on the order of 
10 3 with I, I' ~ 3 mm, by aiming at I , Iq ~ 3 — 30 /xA, 
Z cpw - 50 CI, v p ~ 0.4c, and C c , C' c ~ 2 - 10 fF. 

The resonators were characterized in a 3 He- 4 He dilu- 
tion refrigerator at the base temperature, T < 0.05 K. 
For the characterization, we used a vector network an- 
alyzer choosing an intermediate frequency (IF) between 
100 Hz and 40 kHz. The transmission coefficient S21 of 
A/ 2 resonators were measured by connecting Ports 1 and 
2 in Fig. QJa) to the network analyzer. In the input mi- 
crowave line connected to Port 1, 20-dB attenuators were 
inserted at T = 1 K and at the base temperature. The 
output line had an isolator at T = 1 K, and at T = 4 K, 
an cryogenic amplifier, whose gain at 10-11 GHz was 34- 
38 dB or 40 dB. For measuring the reflection coefficient T 
of A/4 resonators, we employed a directional coupler and 
extra isolators at the base temperature as in Fig. 1(a) 
of Ref. Il2j. Magnetic fields were applied in the direction 
perpendicular to the substrate by means of a supercon- 
ducting solenoid. 

The gains of amplifiers and the insertion losses of at- 
tenuators and cables had been characterized at relevant 
frequencies by independent measurements. The uncer- 
tainties of Pi n , IS21I, an d |r| in this paper are less than 
1-2 dB. 



IV. RESULTS 

Basic parameters of the resonators are listed in Tablefl] 
We obtained fo ~ 10 GHz and Qi on the order of 10 3 as 
we designed. Hystereses were observed in all resonators; 
the critical incident power P c in Table UJ was determined 
experimentally. The values of P c [= (V c /2) 2 /(2Z )} for 
Series A correspond to Io, Iq = 5 — 21 /iA, or the critical 
current density J c = 1 — 2 kA/cm 2 , within the theory^ 
and the mappin g 10 ! 11 in Sec. [Til The estimated critical 
current is on the right order of magnitude, and consistent 
with the dc measurements on test junctions fabricated on 
the same wafer, J c ~ 1 kA/cm 2 . We also estimated the 
critical current by the simulation in Sec. |TTJ taking into 
account the junction capacitance, and confirmed that the 
estimates did not change significantly. The situation is 
similar in Series B. The estimated J c = 0.5 — 0.6 kA/cm 2 
is on the right order of magnitude, and consistent with 
the parameters of the fabrication process, which targeted 
at J c — 0.4 kA/cm 2 . Below, we focus on Resonators Al 
and Bl because within each series, the results were qual- 
itatively similar. 

A. Resonator with a single Josephson junction 

The Pi n dependence of Resonator Al is shown in Fig. 2] 
At different values of P n , the frequency / was swept up 
and down while S21 was recorded. At P n = —101 dBm 
< P c , that is, in the linear regime, the phase of 521 vs. / 
shows a usual rotation [the top right curve in Fig.fJJa)], 
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FIG. 4: (Color online) Dependence on the incident power Pi n 
of Resonator Al. In both (a) and (b), from right to left, 
P in = -101, -93, -85, and -75 dBm. (a) Phase of the 
transmission coefficient S21 vs. frequency. The frequency is 
swept up and down for all curves. For the Pi n = —75 dBm 
curve, the sweep directions are indicated by the arrows. The 
origin of the vertical axis is offset for each curve for clarity, 
(b) Imaginary part of S21 vs. real part of 521. The diameters 
d of the dotted circles range from 0.22 to 0.67 in steps of 0.09. 
The centers are located at ((2/2, 0). 



and the polar plot of Im(52i) vs. Re(52i) falls on the 
largest dotted circle in Fig.[4fb). The center of the circle 
is at (d/2, 0), where d (= 0.67) is the diameter. Note 
that d is equal to the maximum |52i|. It is also related 
to the quality factors by 

Qe/Qu = d~ X - 1. (19) 

Thus, Q e /Q u = 0.5 for Resonator Al in the linear 
regime. The other curves in Fig. 2] are for Pj n > P c , and 
hysteresis is seen as the theory predicts. An interesting 
feature in Fig. [4jb) is that the three curves for P„ > P c 
do not always stay on the largest circle. They move to- 
wards smaller circles as IS21I increases, and the trend is 
stronger at larger P in . This p n dependent quality fac- 
tors are a key to understand Fig. [5) where the measured 
bistable region in grey is compared with Eq. (|6|). 

The grey area in Fig. [5] is determined from many two- 
directional frequency sweeps like Fig. Ufa). If the phase 
difference of S21 between up and down frequency sweeps 
at (/, P in ) was more than 0.1-7T, (/, P in ) in Fig. [5]is grey. 
The pairs of curves in Fig. [5] are Eq. © with different 
values of Qi. The pair with the largest u c is for Qi — Qi 
in the linear regime and associated with the largest circle 
in Fig. 0Jb). The other pairs in Fig. [5] correspond to the 
smaller circles in Fig.@|b) when we neglect the P n depen- 
dence of Q e , which is expected to be much smaller than 
that of d in Resonator Al. The ratio of V c (Qi)/V c (Qio) in 
Fig.[5]is consistent with Eqs. (JTJ) and ([5]), that is, the ratio 
is equal to (Q;/Q ;o )" 3/2 (l - V3/Qi)/(1 - VS/Qiv). We 
treated V c {Qio) as an adjustable parameter, and chose 
its value so that the experimental bistable region stayed 
between the theoretical boundaries. 

The theoretical size of the bistable region in Fig. [5] 
shrinks rapidly with decreasing Qi, and as we have seen, 
Fig. HJb) suggests that Qi of the resonator decreases 
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FIG. 5: (Color online) Bistable region (grey) of Resonator Al. 
The pairs of curves are Eq. © with different loaded qual- 
ity factors Qi. As the power P[ n is increased, the theoretical 
curves with a smaller Qi reproduce the experimental bistable 
region. 



with increasing P; n especially around the lower bound- 
ary, where IS21I takes the maximum value. Thus, p n 
dependent Qi should mostly explain the fact that the ex- 
perimental bistable region is considerably smaller than 
the theoretical prediction for Qi — Qiq. In general, the 
size of experimentally determined bistable region may 
depend on the sweep speed due to finite lifetimes of the 
bistable states. This effect, however, seems to be not sig- 
nificant in Fig. [5]judging from the experiment described 
in the following paragraphs. 

We swept P n instead of / in two directions with 
IF = 40 kHz, which is much faster than IF = 200 Hz in 
Fig. El We chose / = 9.848 GHz and 9.853 GHz, where 
the widths of the bistable region in Fig. [5] are 8.1 dB 
and 5.2 dB, respectively. At both frequencies, the values 
of P; n at which the phase of 52 1 switched [see also the 
arrows in Fig. [3fb)] agreed with Fig. [5] within 0.1 dB, 
despite the large difference in the sweep speed. This is 
the main reason why we presume that the dependence on 
the sweep speed is insignificant in Fig. [5] 

At / = 9.853 GHz, the measurements were repeated 
1.8 x 10 4 times, and the histograms of the switching 
events are shown in Figs.^a) and^b). The distribution 
widths of the switching events are much smaller than the 
width of the bistable region. The histograms can be con- 
verted to the lifetimes r of the bistable states by assum- 
ing that the probability of remaining in the same state 
decays as exp(— At/r), where At is the time and in our 
case, equal to the inverse of IF. This type of conversion 
has been done many times for the conventional switch- 
ing from the zero-voltage state to the voltage state in 
JJs, and the details of the conversion are found, for ex- 
ample, in Ref. 1 2 21 The lifetimes obtained from Figs. [Ha) 
and[(5Jb) are shown in Figs. [Sfc) and^d), respectively. 
Below, we analyze r within a simple model. 

Let us suppose that 

r" 1 cx exp(-P/fc B P), (20) 
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FIG. 6: (Color online) Switching between the bistable states 
at / = 9.853 GHz in Resonator Al. (a) and (b) Histograms of 
the switching events. The horizontal axis is the incident power 
Pin- (c) and (d) Inverse of r vs. P[ n , where r is the lifetime 
of the bistable state calculated from (a) and (b), respectively. 
The curves are the least-squares fittings of Eqs. (|20[) and (|21|l . 



and that as in Ref. |23j, E is associated with the noise. 
When we further assume for simplicity that E is propor- 
tional to the minimum V N , where Vn is the noise voltage 
satisfying Vat + Ue J0JV = V&> V is the source voltage, On 
is the phase noise, and V& corresponds to one of the theo- 
retical boundaries of the bistable region [see the markers 
in Fig. EDJa)], the power dependence of E is written as 



E oc 



'Pr, ~ 



(21) 



because V and V& are in phase when V N takes the mini- 
mum value. Intuitively speaking, the switching between 
the bistable states would occur at Pf, if an experiment 
with no noise were possible. The last assumption is actu- 
ally consistent with the approach based on the Poincare 
section in Fig. 2 of Ref. |2| because the height of the sad- 
dle point in the Poincare section is proportional to the 
minimum |Vjv|. 

The curves in Figs. [He) and[5Jd) are the least-squares 
fittings of Eqs. IJ3DJ and (JUJ). The fittings yielded P b = 
Pi — 0.5 dB and Ph + 0.1 dB. Because P& extends the 
bistable region by 0.6 dB only in total, the above analysis 
of t at f — 9.853 GHz supports our presumption that the 
dependence on the sweep speed is insignificant in Fig. [5] 



B. Resonator with a dc SQUID 

The SQUID modulation of Resonator Bl in the linear 
regime measured at Pj n = —113 dBm is shown in Fig. [7J 
The external dc magnetic flux $ indeed changes Iq, and 
thereby, Lj and fo, periodically. Around &/&o = ±0.5, 
the response was a little hysteretic and fo depended on 
the direction of the flux sweep. In Fig. the flux was 
swept in one direction only, and that is why there is small 
asymmetry at $/<I>o ~ ±0.5. From the amplitude of the 




FIG. 7: (Color online) Resonant frequency fo in the linear 
regime vs. external dc magnetic flux $ normalized by the 
superconducting flux quantum $o = h/(2e) for Resonator Bl. 
The direction of the flux sweep was from left to right. 
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FIG. 8: (Color online) Bistable regions (grey) of Resonator Bl 
when the resonant frequencies fo in the linear regime were 
tuned to 11.16, 11.23, 11.26, and 11.30 GHz by dc magnetic 
flux. The pairs of curves are Eq. (0, where the loaded quality 
factors from left to right are Qi/10 3 = 1.2, 1.2, 1.2, and 1.3. 



fo modulation, we estimate based on the simulation in 
Sec. [II] that at <!>/<I>o ~ ±0.5, the critical current of the 
SQUID becomes about 15% of its maximum value, and 
thus, the SQUID has a normalized loop inductance of 
Lsq/oo/'&o ~ 0.1, where loo is the critical current per 
junction at $ = 0. The above value of Psq^oo/'&o is con- 
sistent with the design. 

For different values of fo, we measured the bistable 
regions. Some of the results are shown in Fig. [5] The 
rightmost data set was obtained when fo was close to 
its maximum value. Regarding quality factors, we found 
two major differences compared to Resonator Al. One 
is that Qe/Qu (~ 5) in the linear regime is much larger. 
This should be due to the layer of sputtered SiC>2 that 
exists only in Series B as we mentioned in Sec. IIIII The 
other difference is that judging from polar plots (data 
not shown) like Fig. 0Jb), Qi does not change very much 
when Pj n is increased from P c . Actually, we estimate at 
-Pin ~ Pc, a larger value of Qi = 1.3 x 10 3 than that in the 
linear regime. The two curves for the rightmost data set 
in Fig.[5]are Eq. © calculated with Qi = 1.3 x 10 3 . They 
reproduce the experimental boundaries of the bistable 
region. 

The other data sets in Fig. [5] for smaller values of fo 
are also compared with Eq. ([S]), where we employed for 



7 




10.5 11.0 
f (GHz) 

FIG. 9: (Color online) Logarithm of the loaded quality factor 
Qi as a function of the resonant frequency fo in the linear 
regime for Resonator Bl. The solid and dotted curves are the 
theoretical predictions based on the simulation in Sec. [IT] The 
solid curve is for R^p — 0, that is, the quasiparticle resistance 
iiqp of the SQUID is not taken into account, whereas the 
dotted curves from top to bottom are for R qp = 20 kf2, 10 kfl, 
5 kf2, and 2 kS7, respectively. 



Qi, the estimated values at p n ~ P c again. The de- 
pendence of P c on /o qualitatively agrees with Eq. ([TJ 
because a smaller f means a larger Ljq. As f is de- 
creased, however, it becomes easier to see some discrep- 
ancies between the experimental boundaries and Eq. (|6j). 
In the leftmost data set, for example, polar plots (data 
not shown) like Fig. [JJb) suggest a moderate decrease 
of Qi with increasing P n above P n ~ P c . Qualitatively 
speaking, the behavior of Qi is becoming similar to that 
in Resonator Al. 

When /o is decreased further down to f < 11 GHz, 
we did not find the bistable region at all, which is a bit 
surprising because according to Eq. ([2]), a larger Ljq is 
more advantageous for observing the bistability, in princi- 
ple. We should not forget, however, that I c /h in Eq. ([2|) 
depends on Qi as well, and in fact, Qi at fo < 11 GHz 
decreases with decreasing fo rather sharply even in the 
linear regime as shown in Fig. [9] where the vertical axis 
is the logarithm of Qi, and Qi was calculated from the 
data at P n = — 113 dBm. We will discuss the origin of 
the change in Qi later in the next section. The lifetimes 
r of the bistable states may be also related to the disap- 
pearance of the bistable region at /o < 11 GHz. As fo 
is decreased, the relevant power range becomes smaller. 
Within the simple model of Eqs. (j2"0)) and (JUJ, smaller 
P; n and p, make E smaller, and thereby r shorter. Thus, 
we cannot exclude the possibility that our measurements 
were not fast enough for observing the bistability at 
fo < 11 GHz. 



V. DISCUSSION 

In the preceding section, we have seen that how the 
quality factor depends on P; n or $ is relevant to the ex- 
perimental bistable region. In this section, let us discuss 



the origin of the dependence. 

Experimental determination of all the quality factors 
in Eq. ([3]) is straightforward in the linear regime. For A/2 
resonators, Qi is obtained from IS21I vs. /, and Q e /Qu 
from Eq. (Til))) by substituting the maximum value of IS21 1 
for d. In Fig. [S] for example, the sharp decrease of Qi at 
/o < 11 GHz is mostly due to the change in Q u . We 
have also confirmed that the frequency dependence of 
Q e is consistent with the theoretical prediction that Q e 
is determined mainly by ljqC c . The unloaded quality 
factor Q u is a measure of the internal loss. In our res- 
onators, both CPWs and JJs can be the source of loss. 
Resonators Al and Bl have completely different Q u in 
the linear regime: Q u — 1.6 x 10 4 and 1.3 x 10 3 , re- 
spectively. As the most probable reason, we have al- 
ready mentioned the difference in the quality of the di- 
electric for the CPWs. The value of Q u = 1.6 X 10 4 for 
Resonators Al is actually still considerably smaller than 
those we obtained for usual CPW resonators without JJ 
fabricated at the same time as Resonators Al. Thus, 
there should be a contribution from JJ as well. 

In the bistable region, it is still possible to estimate the 
quality factors. We have seen that Fig. dfb) suggests a 
strong p n dependence of Q u in Resonator Al. The P m 
dependence is likely to be due to a finite P qp of J J. The 
main reason is that the usual CPW resonators without 
JJ looked almost P n independent in the relevant power 
range. In general, Q u of CPW resonators is known to 
increase with increasing p n at least in the weak-power 
range, where the average energy stored in the resonator 
at u>q is comparable to fru>o or smaller^ This behavior 
was also observed in our resonators in the linear regime. 
Both in Resonators Al and Bl, Q u gained several per- 
cents when P; n was increased by 5 dB from P c — 15 dBm. 
Figure 0Jb) is explained if P qp is a function of Vj and its 
value decreases with increasing Vj, which is quite likely 
(see, for example, Sec. 6.3 of Ref. 24). Although we are 
interested in the case that the JJ stays in the supercon- 
ducting branch with zero dc voltage, it does not nec- 
essarily mean that the ac voltage is also zero. The ac 
voltage is given by Eq. (|11[) . and when it is nonzero, an 
ac current flows through P qp causing a loss. This picture 
is supported by the following quantitative consideration: 
according to the simulation in Sec. [H] the reduction of 
Qi becomes considerable when i? qp is decreased below 
10 3 f2 in Resonator Al, whereas the dc measurements on 
the test junctions suggest that P qp in Resonator Al can 
be as low as 3 x 10 2 fl. It is also consistent with the 
results on Resonator Bl. Because Bl has a smaller J c 
and a much larger superconducting gap, the tunnel bar- 
rier of the junctions is thicker. Thus, P qp must be much 
larger, which explains the fact that p n dependence is 
much smaller. The <3 ) dependence is also explained be- 
cause $ modulates the ratio of P qp / (u)qL jo) . As $ ap- 
proaches ±0.5<I>o, ±1.5$0; ' ' ' 7 Ljo increases, but there 
is no obvious mechanism that $ changes P qp . We have 
also confirmed by the simulation in Sec.[TT]that Qi indeed 
depends on the ratio of R qp /(uJoLjo), and that the ex- 
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perimental Qi vs. /o in Fig.[9]is qualitatively reproduced 
with R qp on the order of 10 3 - 10 4 Q. In Fig. the pre- 
dictions based on the simulation for different values of 
i? qp are also shown by the solid and dotted curves. 

From the viewpoint of the qubit readout and the qubit 
coherence time, the dissipation due to the loss at R qp can 
be unfavorable depending on the circuit configuration, 
and worse than that at CPWs when R qp is closer to the 
qubit. Fortunately, now we know from the experiment in 
Sec. HVl how large R qp has to be in order for the loss at 
i? qp to be negligible. For our resonators in Table HI the 
ratio of R qp / '{ujqLj ) needs to be larger than 10 3 — I0 4 
according to the simulation in Sec. [TXJ This condition 
should be easily satisfied even with Al JJs by employing 
a sufficiently small critical-current density J c . 



region of the resonators were experimentally determined, 
and compared with the theory and simulations. We 
found that the variation of the unloaded quality factor as 
a function of relevant quantities such as the drive power 
and the external magnetic flux, was important for under- 
standing the experimental results. The unloaded quality 
factor is a measure of the internal loss, and the origin of 
its variation was also discussed. 
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